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Abstract: We construct approximate inflationary solutions rolling away from the unstable 
maximum of p-adic string theory, a nonlocal theory with derivatives of all orders. Novel 
features include the existence of slow-roll solutions even when the slow-roll parameters, as 
usually defined, are much greater than unity, as well as the need for the Hubble parameter 
to exceed the string mass scale m s . We show that the theory can be compatible with 
CMB observations if g s / y/p ~ 10~ 7 , where g s is the string coupling, and if m s < 10~ 6 M p . 
A red-tilted spectrum is predicted, and the scalar-to-tensor ratio is bounded from above 
as r < 0.006. The p-adic theory is shown to have identical inflationary predictions to a 
local theory with superPlanckian parameter values, but with the advantage that the p-adic 
theory is ultraviolet complete. 
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1. Introduction 

Many string theorists and cosmologists have turned their attention to building and testing 
stringy models of inflation in recent years. The goals have been to find natural realizations 
of inflation within string theory, and novel features which would help to distinguish the 
string-based models from their more conventional field theory counterparts. The more pop- 
ular categories include brane-antibrane D3/D7 ||, modular ||, DBI § and tachyon- 
driven || inflation (see || for a review). 

In most examples to date, string theory has been used to derive an effective 4D field 
theory operating at energies below the string scale. Since string theory provides a complete 
description of dynamics also at higher energies, it may be interesting to consider a model 



- 1 - 



which takes advantage of this distinctive feature. This is usually daunting since the field 
theory description should be supplemented by an infinite number of higher dimensional 
operators at energies above the string scale, whose detailed form is not known. In the 
present work, we propose to take a small step in the direction of overcoming this barrier, 
by considering a simplified model of string theory invented in 1987 |7j , in which the world- 
sheet coordinates of the string are restricted to the field of p-adic numbers. Scattering 
amplitudes of open string theory can be related to those of the p-adic strings. A great 
advantage in p-adic string theory is that it is possible to compute all amplitudes of its lowest 
state and to determine a simple field-theoretic Lagrangian which exactly reproduces them. 
The result is a nonlocal field theory which is nevertheless sensible in the far ultraviolet. 

The p-adic string resembles the bosonic string in that its ground state is a tachyon, 
whose unstable maximum presumably indicates the presence of a decaying brane, analogous 
to the unstable D25-brane of the open bosonic string theory ||]. Similarly to the bosonic 
string, the potential is asymmetric around the maximum, with one direction leading to a 
zero-energy vacuum, while in the other direction the potential is unbounded from below. 
We will consider whether it is possible to get successful inflation from rolling toward the 
bounded direction. This has been tried before in the context of the open string tachyon, 
and is difficult || because the tachyon potential is not flat enough to give a significant 
period of inflation, and there are no parameters within the theory which can tune the 
potential to be more flat. In contrast, we will show that the p-adic string tachyon can 
roll slowly enough to give many e-foldings of inflation. There are two distinct regions of 
parameter space which allow for successful inflation. There is a region with p = 0(1) for 
which the p-adic field potential is flat and slow roll inflation proceeds in the usual manner. 
However, there is also a region of parameter space with p > 1 for which the potential 
is extremely steep (\rj\ = Mp\V"/V\ may be as large as 10 11 ) but the p-adic scalar field 
nevertheless rolls slowly. This remarkable behaviour relies on the nonlocal nature of the 
theory: the effect of the higher derivative terms in the action is to slow down the field 
sufficiently, despite its steep potential. This new effect manifests itself only in the regime 
where the higher derivative interactions cannot be ignored. It is also interesting that in 
our model the kinetic energy is responsible for driving inflation for a significant number of 
e-foldings, unlike in conventional models of inflation. 

One may worry about the presence of these higher derivative terms, because they are 
usually known to introduce ghosts 1 into the theory. In fact, it is easy to check that for 

1 One may also worry about classical instabilities which usually plague higher derivative theories; gener- 
ically they go by the name of Ostrogradski instabilities (see [pi for a review). These are the classical 
manifestations of having ghosts in the theory: they can have arbitrarily large negative energy, which leads 
to classical instability. Since the nonlocal theory under consideration does not contain any ghosts we also 
do not expect to find such instabilities. One way to see how such theories may avoid the Ostrogradski insta- 
bility argument, valid for finite higher derivative theories, is by noting that one cannot construct the usual 
Ostrogradski Hamiltonian because there is no highest derivative in such nonlocal actions. Also, in arriving 
at the Ostrogradski Hamiltonian, one assumes that all the derivatives of the field (except the maximal one) 
are independent canonical variables. This is no longer true for theories with derivatives of infinite order. 
For instance, it is not possible to independently choose an initial condition with arbitrarily specified values 
of all the derivatives of the field. See ref. fill for a discussion of this point. 
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a scalar field theory if one introduces only a finite number of higher derivatives, then the 
model invariably contains ghost degrees of freedom. The reason why the p-adic string 
action can evade this problem is because it is intrinsically nonperturbative in nature, the 
propagator being modified in such a way as to not contain any poles. In other words 
there are no physical states, ghosts or otherwise, around the true vacuum. This novel 
way of curing the problem of ghosts in higher derivative theories, while retaining some 
nice properties such as improved UV behaviour, was already pointed out in [12] in the 



context of gravity. It was also pointed out in |12| that such theories can exibit interesting 
new cosmological features. For instance, one can obtain nonsingular bouncing solutions by 
making gravity weak at short distances. More recently, such models have also been shown to 



possess inflationary solutions |fl3 |. However, the models of |T^, [13|] are phenomenological, 
while the p-adic action that we consider is an actual (albeit exotic) string theory and 
reproduces many nontrivial features of conventional string theories. 

We start by reviewing the salient features of p-adic string theory in section |2[ In 
section || we show that this theory does not give inflation if the higher-derivative terms in 
the action are ignored. However, near the top of the potential, the energy can be large 
enough to justify keeping all higher derivative terms. In section || we show how to resum 
their contributions and we construct approximate inflationary solutions valid near the top 
of the potential, by solving the coupled equations for the tachyon and scale factor of the 
universe. We give two different approximate methods in this section. In section |5] we solve 
for the fluctuations around this background to determine the power spectrum of scalar and 
tensor perturbations which can be probed by the cosmic microwave background (CMB). 
There we show that it is possible to choose parameters which are compatible with the 
measured amplitude and spectral index, and that the scalar-to-tensor ratio is bounded 
from above as r < 0.006 in this model. We also argue that it can be natural to have initial 
conditions compatible with inflation in the p-adic theory. We give conclusions in section |6|. 
Appendix [A] gives details about the p-adic stress tensor and the approximate inflationary 
solution of the Friedmann equation. Appendix ^ gives mathematical details about the 
incomplete cylindrical functions of the Sonine-Schaefli form. Appendix [C] explains a formal 
equivalence between the dynamics of the p-adic tachyon and those of a local field theory 
with a super-Planckian vacuum expectation value (VEV). 

2. Review of p-adic string theory 

The action of p-adic string theory is given by 

= < I ' d * x (-he-^f + J—^i) (2.1) 
9p J V 2 V + 1 

where □ = — d\ + V 2 in flat space and we have defined 



2 9 2m 2 . . 

and ml = -—^ (2.2) 

" In n 



1 _ 1 p 2 _ 2m; 

— k = — k ~ and fn v = - — 
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The dimensionless scalar field 4>(x) describes the open string tachyon, m s is the string mass 
scale and g s is the open string coupling constant. Though the action ( |2,1[ ) was originally 
derived for p a prime number, it appears that it can be continued to any postive integer 
and even makes sense in the limit p — > 1 [Q. Setting □ = in the action, the resulting 
potential takes the form V = (mf /^p)(^0 2 — Its shape is shown in figure [|. 
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Figure 1: The potential of the p-adic tachyon for p = 19. 

The action ( p.l| ) is a simplified model of the bosonic string which only qualitatively 
reproduces some aspects of a more realistic theory. That being said, there are several 
nontrivial similarities between p-adic string theory and the full string theory. For example, 
near the true vacuum of the theory <p = the field naively has no particle-like excitations 
since its mass squared goes to infinity. 2 This is the p-adic version of the statement that 
there are no open string excitations at the tachyon vacuum. A second similarity is the 
existence of lump-like soliton solutions representing p-adic D-branes [JDJ. The theory of 
small fluctuations about these lump solutions has a spectrum of equally spaced masses 



squared for the modes [15|,[16|, as in the case of normal bosonic string theory. 

It should also be noted that the connection between ( |2.1| ) and the DBI-type tachyon 
actions, which have been widely studied in the literature in the context of tachyon matter 



1 17], is not entirely clear (see for a discussion of the relation between p-adic and 
ordinary strings). In the case of tachyon matter, solutions which roll towards the vacuum 
T — > oo have late time asymptotics T ~ t and hence the tachyon never reaches this point 
|19|| , whereas in the case of the p-adic string the vacuum is at a finite point in the field 
configuration space and homogeneous solutions rolling towards the vacuum typically pass 
this point without difficulty |llj (at least in flat space). In fact, the numerical studies 
of |il~|| found no homogeneous solutions which appeared to correspond to tachyon matter 



(vanishing pressure at late times) . This issue has been considered in the case of cubic string 
field theory in [20]. Related rolling tachyon solutions in cubic string field theory have been 
discussed in 



2 Reference jTl| found anharmonic oscillations around the vacuum by numerically solving the full non- 
linear equation of motion. However, these solutions do not correspond to conventional physical states. 
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It is worth pointing out that one obtains very similar actions to (|2.l| ) with exponential 
kinetic operators (and usually assumed to have a cubic or quartic potential) when quantiz- 
ing strings on random lattice ]23] . These field theories are also known to reproduce several 
features, such as the Regge behaviour J23|] . of their stringy duals. Although our analysis 
focuses on the specific p-adic action, it can easily be applied to such theories as well. The 
connection between p-adic string theory and ordinary string theory on a discrete lattice 
was explored in [18]. 



The field equation that results from ( |2.1| ) is 

e - D l m l <ft> = <ftP (2.3) 

We are interested in perturbing around the solution 0=1, which is a critical point of the 
potential, representing the unstable tachyonic maximum. For odd p one also has another 
unstable point at (ft = — 1, but we will restrict our attention to solutions that start to evolve 
from the omnipresent (ft = 1 maximum. In passing we also note that there is also the stable 
vacuum of the tachyon, at (ft = 0. For both even and odd p the potential is unbounded 
from below. 

It is also worth commenting on the physical interpretation of the fact that the potential 
is unbounded from below. The instability associated with the decay of the "closed string 
vacuum" (ft = to the "true vacuum" (ft = oo is thought to be associated with the closed 
string tachyon instability [||. 

Notice that in the limit p — > 1 the equation of motion ( [2.3| ) becomes a local equation 

U(ft = 2(ftln(ft 

Therefore p — > 1 is the limit of local field theory. In the present work we will also consider a 
very different limit, p S> 1, in which the nonlocal structure of (^]l]) is playing an important 
role in the dynamics. 



3. Absence of Naive Slow Roll Dynamics 



One may wonder whether the field theory (2A) naively allows for slow roll inflation in 
the conventional sense. Naively one might expect that for a slowly rolling field the higher 



powers of □ in the kinetic term are irrelevant and one may approximate (2A) by a local 
field theory. The action (2.1) can be rewritten as 



S 



d A x 



(3.1) 



where we have defined the field x as 



x = xo<ft 



Xo 



p 



hip 



9s V 2(p - 1) 



m. 



(3.2) 
(3.3) 
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and the potential is 



v(x) = ^x 2 

mp 



p 



9s P 2 



1 



XoJ 



P+i 



(3.4) 



In the • • • denotes terms with higher powers of □. The truncation ( |3.1| ) is quite 

analogous to what is usually performed in the literature when one studies inflation from 
the string theory tachyon ||. In the case of the usual string theory tachyon corrections 
involving D 2 and higher are expected, but the full infinite series of higher derivative terms 
is not known explicitly (see, however, |25f | for a calculation of the string theory tachyon 
action up to order <9 6 ). Thus, the standard approach is to simply neglect such terms. We 
will show, however, that under certain circumstances the higher derivative corrections may 
play an extremely important role in the inflationary dynamics. 



Working in the context of the action (3.1) let us consider the slow roll parameters 
describing the flatness of the potential (3.4) about the unstable maximum x = Xo- 
straightforward to show that 



It is 



M 2 



1 



(dV( X )\ 



2 V{xo? \ d X J 
1 d 2 V( X ) 



P V( X o) d X 2 



x=xo 



X=XO 



2 



1 (Mp 

Inp p 2 \ 



m. 



(3.5) 
(3.6) 



Thus one would naively expect that inflation is only possible in the theory (24.) for 
g s M p /(^lnpm s ) <C 1. For p = 0(1) we will see that this expectation is correct, how- 
ever, for p> 1 this intuition is incorrect and that successful inflation can occur even with 
M 2 V~ 1 d 2 V/dx 2 ~ — 10 16 ! The reason for this surprising result is that the dynamics of 
the p-adic tachyon field is set by the mass scale which appears in the kinetic term, m s , 
rather than the mass scale which is naively implied by the potential 



d 2 v 



dx 2 



x=xo 



2(p ~ 1 ) 
lnp 



-m Q 



(3.7) 



Clearly for p S> 1 we have | m 2 | 3> m? a 



4. Approximate Solutions: Analytical techniques 

In this section we construct the approximate solutions for the scalar field and the quasi-de 
Sitter expansion of the universe, in which <j) starts near the unstable maximum (</> = 1) 
of its potential and rolls slowly toward the minimum (0 = 0). Cosmological solutions of 
similar nonlocal theories have also been considered in |24j| . 

We use two different formalisms to construct inflationary solutions. We first devise 



a perturbative expansion in e xt similar to what was carried out in |11] to study rolling 
solutions in flat spacetime. Our second formalism is the analogue of the usual slow roll 
approximation: we assume that the friction term in the □ operator dominates over the 
acceleration term and also neglect the time variation of H. 
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We first discuss the perturbative expansion in powers of e xt . Our starting point is the 
ansatz 

oo 

<Kt) = 1 - J> r e rAt (4.1) 

r=l 

and 

oo 

H{t) =H ~Y^ H r e rXt (4.2) 

r=l 

We have chosen the parameterisation such that at t — * — oo, <p starts from the top of the 
hill where the universe is undergoing de Sitter expansion with Hubble constant Hq. As 
t — ► — oo, e A< — > and all the correction terms, which quantify the departure from the pure 
de Sitter phase, vanish. As t increases, the field rolls toward the true vacuum <fi = 0, in 
fact reaching it at a finite time. Classically, the model admits infinitely many e-foldings 
of inflation, although only the last 60 e-foldings before the end of inflation are relevant 
for observation. This idealized behaviour is an artifact of neglecting quantum fluctuations; 
quantum mechanically the field cannot sit at (j) = 1 for an infinite amount of time. We 
will return to this issue later and show that the inclusion of quantum fluctuations does not 
spoil inflation, which it would in ordinary local field theory if the r] parameter is large. 
Given the ansatz (4.1,4~2|), we can expand the field equation for the p-adic scalar (^ 



and the Friedmann equation as a series in e xt and then determine the coefficients {(j) r , H r } 
systematically, order by order. We will show that this can be done consistently. The zeroth 
order Klein-Gordon equation is trivially satisfied, by virtue of the fact that we start from 
a maximum of the potential. 

4.1 p-adic Scalar Field Evolution 

Let us first find an approximate solution for the scalar field equation of motion Q2,3| ). We 
note that to compute quantities such as d n e xt to the order of interest, it is sufficient to 
use a truncation of ( |4.1| ) and fl4.2j ) : 

<j) = 1 - u - <j) 2 u 2 , H = H - H lU - H 2 u 2 (4.3) 

where we have used the freedom to choose the origin of time to set = 1, and for 
convenience we have defined the new variable 

u = e xt (4.4) 

in terms of which the □ operator takes the form 



□ = -A 2 



2 d 2 ( _ 3H\ d 



U du 2 I A J du 



(4.5) 



We wish to compute quantities such a n n e xt up to 0{u 2 ). This can be done recursively. 
Writing U n <p (where n > 1) as 

(_U) n <f ) = A n ,u + B n u 2 + ... (4.6) 
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and applying another □ operator one finds the following recursion relations for the coeffi- 
cients A n and B n : 

A n -\-\ = A n (A 2 + 3H X) (4.7) 



and 



(4.8) 
(4.9) 



B n+1 = B n (4X z + 6H X) - 3H l XA n 
Equation ( [4.7] ) has the solution 

A n = ai (X 2 + 3H X) n 
with a\ = — 1 (from examination of □(/>) while a suitable ansatz for B n is given by 

B n = h(X 2 + 3H X) n + 6 2 (4A 2 + 6H X) n (4.10) 
The coefficients b\ , 62 can be deduced from (fO|) using the initial values 



#o + A 



and 60 



H + X 



(4.11) 



which follow from explicitly computing Dcp and D 2 (/>. Putting everything together we now 
have 



(-□) n = d n fi - (A 2 + 3H X) n u - 



#1 , x 2 



#o + A 



(A + 3H X) n 



H + X 
u 2 + 0(u 3 ) 



2 (4A 2 + 6H X) n 



(4.12) 



which works also for the case n = 0. Using ( p^ ) 

one can resum the contributions coming 
from all the powers of □ in the exponential operator e _D//m p to give 



-□ I ml 



1 - e^u + [(a - 2 ) e^ 2 - cre w ] u 2 + 



where we have introduced 

A 2 + 3HqX 



Ml 



mf 



4A 2 + 6H X , Hi 
H2 = ^ and a = 



m- 



H Q + X 



(4.13) 



(4.14) 



We conjecture that such resummations are possible for higher order terms as well. Notice 
that ( 4.13j ) reduces to <j> in the limit m p — ► 00, as it should. 

To solve the equation of motion for the scalar field, we must equate ( 4.13| ) to the 
right-hand-side of ([Q|): 



[1 — u — <f>2U 2 ] p = 1 — pu — p 



p-1 



u 2 + 



(4.15) 



matching coefficients for each order in u. The zeroth order equation is identically satisfied, 
as promised earlier, while matching at first order gives 



= mi = J\ 2 +3H )/ml 



P 



(4.16) 
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which, using (|2.2j ) can be rewritten in the form 



A 2 + 3H X = 2m 2 



(4.17) 



which is independent of p. Later we will see that Hq 3> m s is necessary for getting inflation, 
so the solution of (4.17) is approximately 



A 



2m% 
3H~n 



(4.18) 



Finally, matching coefficients at second order gives 



p — 1 



(4.19) 



In summary, by solving the scalar field equation of motion fl2.3|) to second order in the 
expansion in powers of u = e xt , we have obtained two relations, ( 4.17| ) and ( 4.19 ). The 
former determines the parameter A while the latter determines 4>2- 



4.2 The Stress Energy Tensor and the Friedmann Equation 

To complete our approximate solution for the classical background, we must solve the 
Friedmann equation 



3M p 2 



(4.20) 



to second order in it. To find the energy density pj,, we turn to the stress energy tensor for 
the p-adic scalar field. A convenient expression for T^ u was derived in |2(| (see also jp7j) 



m s 

i 



H o / dr ( 

m p Jo 



p + l 



AP+ 1 



+ 



1 

m p JO 



dr I Ue ™p, 



(4.21) 



OL , ^ 



™>p9p Jo 



1 / _tO 

dr I d,,e m p 



(1— r)P 



One may verify that the T^ u is symmetric by changing the dummy integration variable 
r — > 1 — r in the last term. For homogeneous 4>{t) the above expression simplifies, and for 
Tqo we find 



-Too 



ml 



i 



p + i 



1 / tD 

H ^ / dr \Ue m v< 

™ P Jo 



(I-t)D 



+ 



dr d t [e m p(j)\ d t [e m i ■ 



(4.22) 



One can evaluate the above expression term by term, keeping up to 0(e 2Xt ) ~ u 2 . The 
final result reads 



Vnr, = ^ [l - nil + e^) - ^"^^ + «(^ - 1) 



25 2 



l-u(l + e w ) 
m*(p - 1) 



p+l 



+ C(« 2 



2^(p+i; 



(4.23) 
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The 0(u) terms cancel out and matching the coefficients in the Friedmann equation gives 
us the simple results 

Hl = ^ P ^TT) (424) 



and 



H x = (4.25) 



for zeroth and first order respectively. 

The 0(u 2 ) contribution to Too is quite complicated (see appendix [A|) but once we use 



(4.25) it simplifies greatly. Matching coefficient at order Oiu ) in the Friedmann equation 



gives 

Because of our sign convention for H r , the fact that H 2 > means that the expansion 
is slowing as (ft rolls from the unstable maximum, as one would expect in a conventional 
inflationary model. 

Using ( [4,18| ), (|4.24 ) and ( 4.26 ) to compute H2/HQ ~ plnp it is clear that the pertur- 



bative expansion in u breaks down once u ~ 

(plnp)" 1 / 2 (recall that H = Hq — H 2 u 2 ). 
Thus one expects that once u ~ p -1 / 2 then inflation ends. We verify this claim using an 
alternative formalism in the next subsection. 

To summarize, we have determined the five parameters A, 4>2, Hq, Hi and H2 which 
appear in the solutions for (ft(t) and H(t) up to 0(u 2 ) through the equations ( 4.17 - 4~Tgj ) , 



(4.24- 42^ ). As a check of our result, we can take Hq <C m s and compare it to the the 



Minkowski background solution that was found in [11]. For Hq <C m s we see from (4.17) 
that 

A - V2m s ( 1 - -^-) (4.27) 



2\/2m t 

the first term corresponding to the known Minkowski result. We can also compute the 
coefficient <ft>2 in this limit from ( 4.19| ), 

1 



2 = ~ — — (4.28) 



2(p 2 + p + 1) 

This too coincides with the coefficient that was determined in |i Lj for p = 2. 
4.3 The Friction-Dominated Approximation 

In the previous subsections we have constructed an approximate solution for the p-adic 
scalar rolling down its potential by performing an expansion of (ft, H in a power series in u. 
Furthermore, we have shown that once u ~ p -1 / 2 then this solution breaks down (because 
the 0(u 2 ) term in H{t) become larger than the zeroth order term, Hq). At this point 
inflation has ended. Because the equations of motion are complicated, we now verify this 
behaviour using an alternative formalism which does not rely on small u. 

The method is the same as the slow-roll approximation in ordinary inflation, which 
assumes that (ft <C H(ft. To justify it within the p-adic theory, we will provisionally assume 
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that 



A 2 < 3H X 



so that the evolution is friction-dominated in the usual sense. The consistency of this 
approximation will be established when we match the theory to the observables from the 
inflationary power spectrum later, in eqs. fl5l3D , dUp . It follows that it is a good ap- 
proximation to take 

-a^3H d t 
Then the p-adic scalar field equation becomes 

e adt (j) = <j? (4.29) 

where we have defined 

a = ^ (4.30) 

Our procedure is to treat Hq as exactly constant, solve for (j>(t) and compute the energy 
density p^. If p^ is approximately constant then this series of approximations is self- 
consistent and the solution is reliable. Once ps begins to deviate significantly from a 
constant value then the solution breaks down and we conclude that inflation has ended. 

We now proceed to solve ( |4.29| ) for constant Hq. To this end we expand cj>(t) in Fourier 
modes as 

i r+oo 

<f>(t) = -= / dke- ikt <t> k (4.31) 
\l In j-oo 

so that 



e adt <t){t) = / dke- ikt <j) k (4.32) 



2tt 
1 



oo 



dke -ik(t+ a )^ k ( 433 ) 

'2-K J-oo 

<f>(t + a) (4.34) 



The p-adic scalar equation of motion ( 4.29| ) then takes the simple form 

<f>(t + q) = cf)(t) p (4.35) 

It is quite remarkable that this equation admits dynamical solutions. It is straightforward 
to check that ([4.35|) has the solution 

0(f) = e ~ u ® (4.36) 

where u(t) = e xt and A = aT 1 lnp = 2mg/(3Ho), as before. One can easily check by acting 
on fl4.36|) with the full operator — □ = <9f + 3Hodt that indeed the friction term dominates 
as long as u £ 1. 

Performing a Taylor expansion of ( ^4.36| ) about u = gives 

1 9 

(h = 1 — u -\- —u H 
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which reproduces our solution in the small-u expansion in the limit that m s <C Hq 3 (see 
equations 4.3 and 5.15| ). We also have 

e - D / m p(f) = e~ up ^l-pu + ^p 2 u 2 

which again reproduces our previous results (which can be verified by inserting the solutions 
for fj,i, fi2, cr and into equation f4.13| ). 

We now proceed to construct the energy density for <p(t) in this approximation. It is 
straightforward to show that 



-rD/m 2 p 



exp (—up T ) 



-(l-T)a/m£0 



exp I — up 



1-T\ 



We write ( |4.22| ) as 



m; 



[Ti + T 2 + T 3 + T 4 ] 



where 



Ti 

n 



-a /ml 



p+l 
l 



— jt / c?t ( De m p i 
m p Jo 



(l-r)D 

e m p « 



ml Jo 



\ I drd f { e m vd)\dt[e 



(l-T)P 



(4.37) 
(4.38) 

(4.39) 
(4.40) 

(4.41) 
(4.42) 



Using the scalar field equation fl2.3| ) and ( 4.36| ) the first two terms, ( 4.39 ) and ( [4.40 ), are 
trivial 



Ti +T 2 



p — 1 



-up—u 



(4.43) 



p+l 

Since this term is proportional to (jj p+l we identify it with the potential energy of the p-adic 
scalar. 

The next simplest term to evaluate is T4 ( 4.42| ) which gives 



AV 

2m 2 a 



p\np I dr exp [— u {p T + p 1 T )] 
Jo 



It is useful to change the variable of integration tow = up T and cast this result in the form 

A 2 , 



2m 2 . 



f up dw 


r 2n 
pu 


\ — exp 




lu U 


UJ 



(4.44) 



The dw integral can be performed exactly (though not in closed form) in terms of special 
functions. Since T4 is subdominant to a contribution coming from T3 we will not investigate 
the behaviour of (4.44) any further. 



3 We will see later that this is the same as taking the spectral index equal to unity n s 
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We now study T3 ( 4.41| ), This can be written in the form 



2m 2 j „ 



up 



+ 



A 2 
2m? 



1 + 



dujuj exp 

3ff 
A 



-co 



pu 

CO 



21 



up 



exp 



-a; 



a; 



(4.45) 



The dominant contribution to T3 is the one proportional to Hq since the evolution is 
friction-dominated. This term is also larger than T4. The leading contribution to T3 + T4 
is then 

pu 

CO 



T3 + T4 



tip 



dco exp 



-a; 



(4.46) 



where we have used the fact that 3Ho\/(2m 2 ) = 1 which follows from ( [4.17 ) when A 2 <C 



3Hq\. Since T3 and T4 contain time derivatives acting on (f> it is natural to identify (|4.46j ) 
with the kinetic energy of the p-adic scalar. 

Let us study the behaviour of the kinetic energy, equation ( f4.46|) , as a function of u. 
The integral in ( 4.46| ) can be performed in terms of the incomplete cylindrical functions of 
the Sonine-Schlaefli form |}2| (see appendix ^ for a review). 

,21 

■ 27rup 1 / 2 S-i(—u, —up; 2iup 1 / 2 ) 



up 



dco exp 



-co 



pu 



We now study the behaviour of this integral in various limits. We assume that p S> 1 
throughout since the previous method of expanding in u is valid in the case where p ~ 1. 
At very early times, u < 1/p, this integral goes to zero as 



up 



du exp 



-uj 



pu 



21 



up for u < 1/p 



For intermediate times, 1/p <§i u ^ 1/p > , it is a good approximation to treat the upper 
and lower limits of integration as and +00 respectively. In this approximation we can 
write the incomplete cylinder function in terms of a Hankel function of order v = —1 (see 
appendix [B| for details) . The small argument asymptotics of H^l gives 



up 



du exp 



pu 

CO 



1 + 2u 2 p 

ln(2V /2 ) for 1/p < u < 1/p 1 



Finally we consider late times, u > 1/p 1 / 2 . It is still reasonable to extend the integral as 
f™ p dco = J °° dco and hence the integral can still be written in terms of H_\ . This time the 
large-argument asymptotics of the Hankel function are appropriate and one has 

,2- 



pup 
J u 



dco exp 



-co 



pu 

CO 



TTUp 



1/2 -2UP 1 / 2 



for u > 1/p 



1/2 



We have verified these asymptotic expressions numerically. 

We can now write the dominant contribution to in the friction-dominated approxi- 
mation, 



in: 



P-1 , 

Zg'i p - 1 \p + 1 



-up—u 



+ 



up 



dco exp 



-co 



pu 

CO 



(4.47) 
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Figure 2: Plot of the energy density of the p-adic scalar as a function of t, for p = 10 5 . This 
figure shows that is approximately constant for u < p -1 / 2 . The individual contributions from 
the potential and kinetic energy are also shown. 

The first term, proportional to e~ up , represents the potential energy and the second term 
represents the kinetic energy. Using our previous analysis of the kinetic energy the be- 
haviour of p& as a function of u (assuming p S> 1) is clear. At early times u < 1/p 
the potential energy dominates and we have p^ = rn^/{2g 2 ). At intermediate times 
1/p <C u < 1/p 1 / 2 the potential energy goes to zero and the kinetic energy dominates 
and we have p = m^/(2g 2 ). At late times u > 1/p 1 / 2 and p^ damps to zero as e ~ 2upl/2 . We 
have verified this behaviour numerically. Figure [2] shows the behaviour of p^ as a function 
of t, verifying that p^ is approximately constant for u < p' 1 / 2 . The time evolution of both 
the potential energy and the kinetic energy contributions to p^ are shown, demonstrating 
that the latter dominates in the interval p^ 1 < u < p^ 1 ! 2 . In this figure we have taken 
p = 10 5 for illustrative purposes. 

We conclude that for u < 1/p 1 / 2 the scalar field energy density p^ is approximately 
constant and inflation proceeds. At u = 1/p 1 / 2 the energy density of the p-adic scalar 
begins to decrease quickly and the analysis of this subsection is no longer applicable. We 
conclude that inflation ends, roughly, when u = 1/p 1 / 2 . 

It is quite interesting that during the intermediate phase p^ 1 < u < p^ 1 / 2 it is in fact 
the kinetic energy which is driving inflation, rather than the potential energy. This is quite 
different from what occurs in a local field theory. 

5. Fluctuations and Inflationary Predictions 

In this section we consider the spectrum of cosmological fluctuations produced during p- 
adic inflation. The full cosmological perturbation theory for the p-adic string model (|2,1|), 
which should include metric perturbations and also take into account the departure of the 
background expansion from pure de Sitter, is complicated and beyond the scope of the 
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present paper. We leave a detailed study of these matters to future investigation 
To simplify the analysis we will neglect scalar metric perturbations as well as deviations 
of the background metric from pure de Sitter space during horizon crossing. In standard 
cosmological perturbation theory these approximations reproduce the more exact results 
to reasonable accuracy and we therefore assume that the situation is similar for the action 
(|2J 



5.1 p-adic Tachyon Fluctuations 

We are approximating the background dynamics as de Sitter which amounts to working in 
the limit u — * so that 

tt2 — tt2 _ m t P - 1 /e -,\ 

(ft = <fto = 1 (5.2) 

We expand the p-adic tachyon field in perturbation theory as 

<t>(t,x) = <ft {0) (t) + 8<ft(t,x) 
= 1 + 5<j)(t, x) 

The perturbed Klein-Gordon equation ( |2.3D takes the form 

e- n l m l8<l) = p5<j) (5.3) 



(Inhomogeneous solutions in p-adic string theory have also been considered in |2lf| .) One 
can construct solutions by taking 8(ft to be an eigenfunction of the □ operator. If we choose 
8<ft to satisfy 

-U86 = +B S(j) (5.4) 



then this is also a solution to (5.3) if 

B = rriplnp = 2m 2 s (5-5) 



where in the second equality we have used (|2.2| ) 



The solutions of (5^4) are well known. However, in order to make contact with the 
usual treatment of cosmological perturbations we need to define a field in terms of which 
the action appears canonical. This presents a serious difficulty because, in general, there 
is no local field redefinition which will bring the kinetic term eft ^1 — e~ D ^ m p^j (ft into the 
canonical form (ftO(ft. (One might imagine simply truncating the expansion in powers of 
□ as we have described in section ||, however, the higher order terms are not negligible in 
general.) Fortunately, for fields which are on-shell (that is, when (|5,4j ) is solved) the field 
obeys 

1 _ e -a/"$) 66= (l- e B l m l 
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1 - e 



B/mi 



(-B) 



p-1 
2mA 



Thus, for on-shell fields the kinetic term in the Lagragian can be written as 

/:o„-shell = 4^(l-e- D /^)0+--- 

9p * 

mi p — 1 f , 

1 

= -¥>□¥> + ••• 
In (15.61) we have defined the "canonical" field 



(5.6) 



99 = A(f> 



(5.7) 



where 



.4 



m s p 



(5.8) 



The field has a canonical kinetic term in the action, at least while Q5.4p is satisfied. 
Notice that </? is distinct from the field x ( see ( |3.2| )) which we introduced in section ^. The 
field x corresponds to the canonical field which one would naively define when neglecting 
terms 0(D 2 ) and higher in the action (as is typical in studies of tachyonic inflation) while 
<p is the appropriate definition of the canonically normalized field when taking into account 
the infinite series of higher derivative corrections. 

Now, let us return to the task of solving (|5.4| ), bearing in mind that 5<p = A5(p is the 
appropriate canonically normalized field. We write the quantum mechanical solution in 
term of annihilation/creation operators as 



5ip(t, x) 



d 3 k 

(2vr) 3 /2 



and the mode functions (pk(t) are given by 



a k ip k (t)e ikx + h.c. 



<Pk(t) 



a 3 H 

where the order of the Hankel functions is 



" e f (H-l/2) H (l) ( k 



aHr 



(5.9) 




(5.10) 



and of course a = e Hot . In the second equality in (5.10) we have used ( |5.5| ) and (2.2). In 
writing ( p.9[ ) we have used the usual Bunch-Davies vacuum normalization so that on small 
scales, k 3> aHo, one has 

a" 



-1 



2k 
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which reproduces the standard Minkowski space fluctuations. This is the usual procedure 
in cosmological perturbation theory. However, we note that the quantization of the theory 
fl2.1| ) is not transparent and it might turn out that the usual prescription is incorrect in 
the present context. We defer this and other subtleties to future investigation. 4 On large 
scales, k <C clHq, the solutions (^S|) behave as 



H ( k \ 



3/2-1/ 



aH J 



V2k 3 

which gives a large-scale power spectrum for the fluctuations 



Pi 



5<p 



H V ( k x n °~ x 



2ir J \aHo 
with spectral index 

n s — 1 = 3 — 2v 

From ( ETIcD it is clear that to get an almost scale-invariant spectrum we require m s <C Hq. 
In this limit we have 

4 (m s \ 2 

"•-"-iUJ (5 - u) 

which gives a red tilt to the spectrum, in agreement with the latest WMAP data [29]. For 



n s = 0.95 one has m s = 0.2Hq. Comparing ( |5.9| ) to the corresponding solution in a local 
field theory we see that the p-adic tachyon field fluctuations evolve as though the mass- 
squared of the field was — 2m 2 s which may be quite different from the mass scale which 
one would infer by truncating the infinite series of derivatives: d 2 V/d\ 2 {x = Xo) ( see ec l- 
( |3.7| )). The fact that Hq > m s is an unusual feature; we will comment on it below. 

It is worth pointing out that we have constructed solutions of a partial differential 
equation with infinitely many derivatives, for which we are free to specify two initial data; 
to obtain eq. |5.9| we fixed these using the Bunch-Davies prescription. Precisely the same 
result was obtained when inhomogeneous solutions were studied in |^]. Partial differential 
equations with infinitely many derivatives constitute a new class of equations in mathe- 
matical physics about which little is presently known. In particular, it is not clear how 
to pose the intial value problem for such equations. We conjecture that the most general 
solutions of ( |2.3| ) are specified by two initial data (for example 0(0, x) and </>(0, x)), just like 
equations containing only one power of □. See |}0| for mathematical work on constructing 
solutions of equations with infinitely many derivatives. 

5.2 Determining Parameters 

We now want to fix the parameters of the model by comparing to the observed features of 
the CMB perturbation spectrum. There are three dimensionless parameters, g s , p and the 
ratio m s /M p . The important question is whether there is a sensible parameter range which 



4 Our prescription for choosing the vacuum has the property that in the local limit p — > 1 the cosmological 
fluctuations are identical to the well-known solutions in local field theory. 
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can account for CMB observations, i.e., the spectral tilt and the amplitude of fluctuations. 
Using ( 4.24D in (5.11), we can relate the tilt to the model parameters via 



\ n ,- ll = *j£pl(^) 2 d _ feY- ^J,, (5,2) 



P 2 \ m s J S \Mp) p 2 \n s — 1 

Thus one can have a small tilt while ensuring that the string scale is smaller than the 
Planck scale, provided that g 2 /p <C 1. Henceforth we will use ( |5.12|) to determine m s /M p 
in terms of p, g s , and \n s — 1| = 0.05. All the dimensionless parameters in our solution, 
m s /Ho, X/m s , 4>2 and H2/m s , are likewise functions of n s — l,p and g s . From ( \i .24] ) and 



( |5.12[ ) we see that for p » 1, 



^ - v^ P ^ - ^ ^ IVW^M (5-13) 

-no m s \ p m s Z 

It may seem strange to have H exceeding m s since that means the energy density ex- 
ceeds the fundamental scale, but this is an inevitable property of the p-adic tachyon at 
its maximum, as shown in eq. ( 4.23j ). This is similar to other attempts to get tachyonic 



or brane-antibrane inflation from string theory, since the false vacuum energy is just the 
brane tension which goes like ra\jg s . 

Next we determine A/m s , where A is the mass scale appearing in the power series in 
e xt which provides the ansatz for the background solutions. Consider eq. ( f4.16| ) for A in 
the Hq 3> m s limit. The positive root for A gives 



A /In, — II , 



From (|5l3D and ( [U4] ) we find that 3H \ » A 2 which means that the evolution is friction- 



dominated in the usual sense. As for fa, from eq. ( pif ) it follows that 



2 = (5-15) 



2p i+P - l 
with 

9* |n s -l|/3 + 0(|n s -l| 2 ) 

so that (3 <C 1- Notice that for p> 1, we have (f>2 — — 0.5p~^. Finally we have H2 which 
from eq. ( [4.26 ) is given by 

El = PP^A ^ ( 5.i 6 ) 

To go further, we must impose the COBE normalization on the amplitude of the density 
perturbations, to show that it is possible to satisfy all the experimental constraints while 
keeping m s /M p < 1. The latter requirement is usually needed for the validity of any 
4D effective description of string theory. In a compactification of the d extra dimensions 
whose volume is of order Vd ~ TnJ , we would have M p ~ m s whereas more generally 
M 2 = m 2+d Vd. The 4D effective theory would normally need to be supplemented by 
higher dimension operators if was small compared to m~ d . 
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5.3 Curvature Perturbation and COBE normalization 



In order to fix the amplitude of the density perturbations we consider the curvature per- 
turbation £. We assume that 

C ~ — rSip 

<P 

as in conventional inflation models. To evaluate the prefactor H/ip we must work beyond 
zeroth order in the small u expansion. We take (f) = 1 — u to evaluate the prefactor, even 
though the perturbation 5\ is computed in the limit that <p = 1. This should reproduce 
the full answer up to 0{u) corrections. The prefactor is 

H H 



AXu 
P 



1 



1 



Jl.s 



1 u 



We should evaluate u at the time of horizon crossing, defined to be approximately 
60 e- foldings before the end of inflation t en d, assuming that the energy scale of inflation 
is high (near the GUT scale). We must therefore estimate t e nd- We have shown in the 
last subsection that inflation ends when u ~ 

l/pi/2 From eqs _ ( p^p^) 

we see that 

Hq/\ = 2/\n s — 1|; therefore we can write the scale factor a(t) = e Hot in the form 



a(t) ^ u(i) 2/l " s-11 



so that a* = e 60 a e nd corresponds to 



(5.17) 



-30|n s 



V 



1/2 



The power spectrum of the curvature perturbation is given by 

t 2 / 7 \ n s -l 



c 



Sip 



A 2 



k 



aH { 



where the amplitude of fluctuations Ar can now be read off as 



2 e 60|n s 



^ 3tt 2 p \n 



-1| 



II 3 



(5.18) 



(5.19) 



(5.20) 



As an example, taking n s 
10~ 10 by choosing 



0.95 one can fix the amplitude of the density perturbations 



A c - 



9s_ 



0.48 x 10" 



(5.21) 



To get a more general idea of how the inflationary observables constrain the parameters 
of the model, we will allow n s to vary away from the value 0.95, which is a fit to the WMAP 
data under the assumption that the tensor contribution to the spectrum is negligible. 
Setting A 2 - = 10" 10 and using (|5.20[) gives and expression for g s in terms of p and \n s — 1 1 



3tt 2 



-30|n s -l| 



n. s 



II 3 / 2 x lO" 5 



(5.22) 
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Combining |[2|) with fTl2|) , we also obtain 




/' + 1 e - 3 °K-l|| ns _ l| x 10~ 5 (5.23) 



We graph the dependence of g s and m s /M p on p for several values of the spectral index in 
figures |3] and |j. We see that the string scale is bounded from above as m s /M p < 0.94 x 10 -6 
and that for typical values of p, n s it is close to m s /M p = 0.61 x 10~ 6 . We also see from 



( 5.22 ) that g s is unconstrained and that g s , p are not independent parameters. If we wish 
to take g s ~ 0(1) then we must choose extremely large values p ~ 10 14 . If we restrict 
ourselves to the perturbative regime g s < 1 then this places an upper bound on p: 

p<4.3x 10 14 
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5.4 Comments on Slow Roll and the Relation to Local Field Theory 

It is quite remarkable that our predictions for inflationary observables and also the solutions 
for x, H are essentially identical to the results from local field theory with potential V(ip) = 
g s (ip 2 — v 2 ) 2 /4 where g s v 2 = 2m 2 . (see appendix [C] for a detailed comparison). Indeed, the 
dynamics of the p-adic tachyon are fixed by the mass scale in the kinetic term, m s , rather 
than by the naive mass scale */ —d 2 V/dx 2 (see section |3j), which may be much larger than 
m s . 

It is an interesting feature of this theory that the canonical p-adic tachyon, (f>, can 
roll slowly despite the fact that, working in a derivative truncatation (as in section |3|), one 
would conclude that the tachyon has an extremely steep potential. To see this we first 
define the Hubble slow roll parameters e H , r/ H by 

e H 

£h - Vh 



2ikf2 H 2 



Hip 



(5.24) 
(5.25) 
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These are the appropriate parameters to describe the rate of time variation of the inflaton 
as compared to the Hubble scale. Using the solution <p = 1 — u (recall that <p = Acf), 
A = m s p/{\f2 g s )) we find that 



IP+l^Oln-ll, 

2 p 1 



II 



??.. s . 



II 



(5.26) 
(5.27) 



We see that the Hubble slow-roll parameters, as defined above, are small. This means 
that the p-adic tachyon field rolls slowly in the conventional sense. One reaches the same 
conclusion if one defines the potential slow roll parameters using the correct canonical field, 
which is tp fl5.7|): 
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}_dV_\ 
V~dp) 

, 1 d 2 V 







p V d^ 2 
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1, 

— n, 
2 i s 



(5.28) 
(5.29) 



On the other hand, consider the potential slow roll parameter which one would naively 
define using the the derivative truncated action ([0) : 
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]_dvy 

Vd^J 
1 d 2 V 







vd x 2 



X=X0 



X=X0 



1\ 1, 



hip 



II 



(5.30) 



(5.31) 



where in (5.31) we have used equations ( |3.6| ) and ( ^.12 ). We see that ( ^.31| ) can be enor- 
mous, though the tachyon field rolls slowly. Taking the largest allowed value of p, p ~ 10 14 , 
and n s S 0.95 we have M 2 V~~ x \d 2 V / d\ 2 \ ~ 10 11 ! 

Since large values of p are required if one wants to obtain g s ~ 1, it follows that it is 
somewhat natural for p-adic inflation to operate in the regime where the higher derivative 
corrections play an important role in the dyanimcs. However, in the regime where p ~ 1 
(corresponding to very small coupling g s ) this novel features is not present. For example, 
with p = 3, n s = 0.95 one has g s ~ 10~ 7 and MlV~ 1 \d 2 V/dx 2 \ ~ 0.05 and so the slow-roll 
dynamics are not surprising. 

5.5 Tensor Modes 



Since the p-adic stress tensor ( 4.22 ) does not contribute any anisotropic stresses up to 
first order in perturbation theory it follows that the first order tensor perturbations of the 
metric do not couple to the first order tachyon perturbation (see, for example, In 
fact, the action for the tensor perturbations is given by 



M 2 
p 



grav 



1 



dtd 6 xa{t) 6 -d^hijd^h 13 
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The standard procedure gives a power spectrum for the gravity waves 

k ^ nT 



P T = A\ 



aH 



with amplitude 



We would like to compare this to the power in scalar modes ( |5.20| ). Defining the tensor- 
to-scalar ratio in the usual way 

J_ A 2 

— 100 T 



2 



±A2 

25 



we find that 

1 <? 



2M* H% CH 



(5.33) 



which reproduces the usual result from local field theory. Using ( |5.26| ), we can evaluate r 
as a function of p and n s — 1 

r= lp+l e -60|n.-l| | j| (534) 

2 p 

It is easy to see that r is maximal when p = 2, n s = 1 — 1/60 and hence it follows that the 
scalar-tensor ratio is bounded from above as r < 0.006, which is very small. 

5.6 Comments on Initial Conditions 

As we have previously noted, our classical solution <p(t) sits at the unstable maximum of 
the potential for an infinite amount of time, thus this model would seem to admit infinitely 
many e-foldings of inflation. Of course, this cannot be the case quantum mechanically and 
one expects quantum fluctuations to displace <fi from the false vacuum and cause it to roll 
down the potential (we have assumed that this rolling takes place towards = 0, rather 
than down the unbounded side of the potential). As a consistency check we note that 

= (5.35, 
0n A 



We should compare this to -u* (see equation 5.18| ) , the distance the field has rolled classically 



at horizon crossing. It is straightforward to show that 

10~ 5 |n s - 11 



where we have used (5.21). Since ((<50) 2 ) 1//2 <C u+ for all parameter values it follows that 



the de Sitter space fluctuations (which are present as u — > 0) will not displace the field far 
enough from the maximum of the potential to have any significant effect on the number of 
observable e-foldings of inflation, although they prevent inflation from being past-eternal. 
We note, however, that if one incorporates thermal fluctuations (and initial momentum) 



- 22 - 



then this model may suffer from problems related to fine tuning the intial conditions as in 
small field inflationary models |33|. However, it is not clear if these objections apply to our 
model for several reasons. The first reason is that the dynamics of this theory is peculiar 
and it is not clear how (or if) the phase space arguments of |33|] apply. The second reason 
is that it is not clear what initial conditions for the field <j> are most natural from a string 
theory perspective. Finally we note that since the field cp rolls a distance A > M p in field 
space, our model is not a "small field" model in the conventional sense. 5 

If p-adic superstrings exist, it might also be possible to justify the initial conditions for 
inflation by having topological inflation |p^l , if the tachyon potential is symmetric about 
the unstable maximum. This distinction exists between the tachyon of the open bosonic 
string H (describing the instability of D25 branes), and the tachyon of unstable branes in 
superstring theory [pip . Any realistic extension of the model should have a potential which 
is bounded from below, and if it is supersymmetric, the minima should be at zero, hence 
the additional minima will be degenerate with the one at 4> = 0. The existence of domains 
of the universe in the different minima ensures that there will be regions in between where 
inflation from the maximum of the potential is taking place, so long as the minima are 
discrete and not connected to each other by a continuous symmetry. 

As we have noted previously the fact that the potential V(4>) is unbounded from below 
is thought to be a reflection of the closed string tachyonic instability of bosonic string 
theory. If this conjecture is correct then the addition of supersymmetry should indeed lead 
to a symmetric potential for the p-adic tachyon which is bounded from below, as we have 
suggested above. 



6. Conclusions 

In this paper we have constructed for the first time approximate solutions of the fully 
nonlocal p-adic string theory coupled to gravity, in which the p-adic tachyon drives a 
sufficiently long period of inflation while rolling away from the maximum of its potential. 
In our solution, the nonlocal nature of the theory played an essential role in obtaining 
slow-roll, since with a conventional kinetic term the potential would have been too steep to 
give inflation. One of the novel features of this construction is that the Hubble parameter 
is larger than the string scale during inflation, a condition which would usually invalidate 
an effective field theory description, but which is consistent in the present context because 
of the ultraviolet-complete nature of the theory. 

We found that the experimental constraints on the amplitude of the spectrum of scalar 
perturbations produced by inflation require a small value of the string coupling g s , and 
can be consistent with a large range of values of the parameter p, 1 < p < 10 14 . The 
regime p 3> 1 is interesting because it exhibits qualitatively different behavior relative to 
conventional inflationary models: slow roll despite the potential being steep, and inflation 
being driven by the kinetic as well as potential energy of the field. This regime is also 



3 The skeptical reader might have reservations about the validity of our analysis since A > M p . We note 



that the action (2.1) is not a low energy effective field theory and hence we believe that we are justified in 



using this action even for super-Planckian symmetry breaking scale. 
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interesting because it corresponds to g s = 0(1) and hence appears more natural from 
a string theory perspective. Since the p-adic string is not construed to be a realistic 
model by itself, it may not be very meaningful to question how natural such values might 
be. However, it may not be unreasonable to think of real strings as being composed of 
constituent p-adic strings because the Veneziano amplitude of the p-adic theory is related 
to that of the full bosonic string by A' 1 = Y\ p A p where the product is over all prime 
numbers p. Thus it may not be unreasonable to expect similar behavior to the large-p 
results from a more realistic model. 

The model predicts a red spectrum, in agreement with the latest WMAP data, whose 
tilt is related to the ratio of the string scale to the Hubble rate during inflation via H/m s = 
2/-v/3|n s — 1|. For n s = 0.95 this gives H/m s = 5. This is in contrast to most stringy 
models of inflation which require H < m s in order for the effective field theory to be valid. 
We find the bound r < 0.006 on the tensor modes. It has been estimated that future 
experiments could eventually have a sensitivity of r ~ 6 x 10 -5 |36|] and hence the tensor 
components may in fact be observable. 

We noted that the p-adic model succeeds with inflation where the real string theory 
tachyon fails. But our analysis makes it clear that this could be due to the failure to keep 
terms with arbitrary numbers of derivatives in the action. The effective tachyon action of 
Sen [O] is a truncation which keeps arbitrary powers of first derivatives but ignores higher 
order derivatives, which were essential for obtaining our solution. Thus the new features 
for inflation which we find in p-adic string theory could also be present in realistic string 
theories, if we knew how to include the whole tower of higher dimensional kinetic terms. 
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Note Added 



Upon completing this paper a related work appeared [38] in which an alternative normal- 
ization for the fluctuations of the p-adic scalar was proposed. Motivated by this work, we 
reconsidered our choice of normalization and concluded that ( |5,7l , |5,8| ) is the most appro- 
priate definition of a "canonical" field. Notice, however, that our field ip differs from the 



definition of a canonical field which was proposed in [38]. In [38] a field redefinition is 
advocated which puts the stress tensor T^ v into canonical form, though this definition does 
not have canonical kinetic term in the action. Though we believe that our definition is 
more natural, we stress that in the case of current interest, p-adic inflation, the distinction 
does not generate any significant quantitative difference because our definition differs from 

/hip (which is less than an order of magnitude for 



that of [38] by a factor proportional to 
the values of p which we consider). We are grateful to J. Lidsey for sending us a draft of 
his manuscript prior to publication and also for interesting and enlightening discussions. 
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A. The Stress Energy Tensor and the Friedmann Equation 

Here we compute the 0(u 2 ) term of the approximate solutions 

<f> = 1 — U — 4>2U 2 

H = H - H 2 u 2 
We write the different terms appearing in the energy density as 



f oo 



m s 



p+l 



+ 



1 



1 / _ tO 



(l--r)P 



+ 



™% Jo 



1 / 

dr dt e m p < 



_(i-t)d 

d t ( e m p . 



in 



2 (Ti + T 2 + T 3 + T 4 ) 

5p 



(A-l) 



where Ti, T 3 and T4 are defined as in ( [4. 39 -4. 42). We now evaluate this expression 
term-by-term. 

We first consider T\ + T 2 , which can be written as 



T l + T 2 = ^-±<t? +1 
p+l 



using (HH). Using we have 

p-1 



Ti + T 2 



+ (-p + l)u + (p - 1) 



2 + | 



u 2 + 0(u 3 



p+l 

To evaluate T3, T4 we note that 

e-ro/m 2 ^ = 1 _ e T^ u _ 2e rH2 u 2 + . . . 

e -(l-r)D/ m ^ = 1 _ e (l-r) Wu _ ^ 2e (l-r)r M2n 2 + . . . 



(see equation 4.13 with a = 0). It is straightforward to show that 

□ e - TD/m ^ = e T/Xl (A 2 + 3# A)u + 2 e TM2 (4A 2 + 3H X)u 2 + ■■■ 
We are now in a position to compute T 3 . Using ( |A-3j ) and ( |A-4| ) we have 
1 Z* 1 

T 3 = —z \ dr \{X 2 + 3H \)e 7 
mi Jo 



+ [-(A 2 + StfoAK 1 + 2 (4A 2 + 6H X)e T ^] u 2 + 
The dr integrals are trivial to perform using the identity 
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dre c 



1 



We find that 



T 3 = (p- l)u + [-pui + (e^ 2 - l)<j> 2 ] u 2 + C(w 3 



(A-2) 



(A-3) 



(A-4) 
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We now consider the integrand of T4: 

(i-Qq 



8 t (e ^<PJ d t (e ~^<PJ X 2 u 2 d u (1 - e^u) d u (l - e (1 ~ T)m u) = AVe w (A-7) 
The dr integral is trivial and gives 
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It is straightforward to sum up the various contributions to Too- We find that 
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The fact that the coefficient of the 0{u) term is zero verifies that H\ = 0. We now solve 
the Friedmann equation 



3H 2 = 



M 2 
p 



Ps 



noting that 



H 2 = Hq — 2H H 2 u 2 + 0(u 3 ) 
Matching the coefficients at order u° gives 



H 2 = 



1 m s p — 1 



6</ 2 M p 2 p+l 
as before. Matching the coefficients at order u 2 gives 
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In the first equality we have used that fact that 

4> 2 e^ = pfo-|(p-l) 

which follows from the second order Klein-Gordon equation and in the second equality we 
have used 

A „H X 

— n + 3 — 7T- = hi = hip 



mi 



mi 



"p "p 

which follows from the first order Klein-Gordon equation. Finally we arive at the result 
for H 2 : 



Ho = 



m s pX 1 p lnp / m 



4 5 2 M 2 m 2 8g 2 s p-l \M P 
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B. Incomplete Cylindrical Functions of the Sonine-Schaefli Form 

The Bessel function can be represented by a contour integral of the Sonine-Schaefli form 



Ju(z) 



1 (Z\> 

2tt7 V2 



c+ioo 



uj v 1 exp 



4oj 



du> 



(B-l) 



as long as Re(z^) > —1. In ( |B-1| ) c is an arbitrary positive constant. The incomplete 



cylindrical function, S u (r, s; z) generalizes (B-l) to arbitrary limits of integration 



It follows that 
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Taking the limits r — > and s — > +oo this integral can be written in terms of Hankel 
functions of imaginary argument, K u (z), as 

-v2 ■ 



to v 1 exp 



o 



-UJ 
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(B-4) 



The function K v {z) is real- valued for real z and is related to the usual Hankel function 
as K u (z) = (iir /2)e t7TU ^ 2 Hjp(iz). Using the well known large-argument asymptotics of the 
Hankel functions one may show that 



ui v 1 exp 
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Auj 
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for z>l. 



C. Comparison to Local Field Theory 

In this appendix we perform a detailed comparison of our results for the action (^4|) to 
the theory 



d 4 x 



-d^d^ + v(i/>) 



(C-l) 



with potential 
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where we have defined m 2 = 2g s v 2 . We are interested in obtaining inflation near the unsta- 
ble maximum ifi = 0. The flatness of the potential is parameterized by the dimensionless 
slow roll parameters 
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Unlike in the p-adic theory we do not need to distinguish between e H ,r] H and ev,r/v (see 



equations 5.24-5.31) and hence we drop the subscripts on the slow roll parameters. For 
ip <C v the e parameter ( |C-3| ) is automatically small while 

(M v 



which is small compared to unity for v S> M p . The fact that the symmetry breaking scale 
is large compared to the Planck scale may be reason to doubt the validity of the field theory 



( C-l ). However, for our purposes this is irrelevant. 

It is instructive to consider solving the equations of motion for this theory using the 
formalism of section f|. We begin by speculating solutions of the form 

0(e 3Xt ) 



1> 



ve xt + tp 2 ve 2Xt + 4) 3 ve 3Xt 



H = Hn 



He 



H 2 e 2Xt + 0{e 



3Xt\ 



(C-5) 
(C-6) 



and solving order by order in e xt . The ansatz ( p-5 > C^6| ) is analogous to ( |4.1|J4.2| ) since in 
both cases the (classical) field spends an infinite amount of time at the unstable maximum, 
driving a past-eternal de Sitter phase, before rolling towards the true minimum of the 
potential. We suppose that e xt <C 1 so that ip <C v initially. 
The Klein-Gordon equation is 



± + 3H± 



Plugging in the ansatz ( p-5Q we obtain 



2m 



+ 2m; 



A + 3i^ A = 2mi 



(C-7) 



at first order in e xt . This result is identical to (4.17). At second and third order respectively 
we find 



V'2 
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The Friedmann equation is 
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At zeroth order in e xt we obtain the familiar result 
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At first and second order we find that 
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In writing the second equality in ( C-10| ) we have used A 2 + 3Hq\ = 2m 



2 



Though a complete treatment of inhomogeneities including metric perturbations and 
nonzero slow roll parameters is straightforward in this context we choose to analyze this 
theory using the same approximations as we used in section |H| in order to make more 
explicit the comparison between the two theories. The perturbed Klein-Gordon equation 
(neglecting metric fluctuations and in the limit e xt — ► 0) is 



d 2 t + 3H d t - - 2m 2 ) ^ = 



where a = e Hot . The large scale solution is 



\5ipk\ 



V2k 3 
where 



H ( k \ 



3/2-1/ 



aHo) 




Near scale-invariance of the spectrum requires m s -C Hq. In this limit we have 



Amj „ ( M, 



2 



^-l = -3^ = -8( — ) (C-ll) 



which is identical to ( 5.11j ). This result also reproduces the full calculation incorporating 
metric perturbations: n s — 1 = 2r\ — 6e = —8Mp/v 2 . 

We can use ( |C- 1 1| ) to write the dimensionless quantities Ho/m s , H2/m s and \/m s in 
terms of \n s — 1|. The solution of ( |C-7| ) can be written as 



m. 



A _ . I\n s — 1| 



which is identical to ( |5.14| ). For this solution we of course have 3XHq S> A 2 so that the 
evolution is friction dominated. It is also straightforward to show that 



o 



m s yj3\n s - 1| 
H2 = Hq 

We see that the solutions H, x f° r the theory ( |C-1| ) are identical to those of the theory 
( |2.1| ) up to order e xt . At order e 2Xt and higher, however, the dynamics of the two theories 
differs. 

From ( |C-3 ) one can check that e = 1 at ip ^ v so that it is a good approximation 



to suppose that inflation ends at e xt = 1. It is straightforward to impose the COBE 
normalization V/(eMp) = 6 x 1CP 7 for this model (which will impose g s <C \n s — 1| 2 ), 
however, this is unnecessary for our purposes. It is clear that the inflationary dynamics 
and predictions predictions of the theory ( p-l[ ) are identical to those of the theory fl2.1|). 
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